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ABSTRACT 

Two  dimensional  distortions  of  a  panel  structure  due  to 
welding  are  analyzed  by  the  finite  element  displacement  method 
with  an  assumption  of  elastic  deformation  during  welding 
process.   Construction  of  stiffness  matrices  relevant  to  the 
welding  problem  is  demonstrated  for  the  one  and  the  two 
dimensional  cases.   Computer  programs  are  presented. 

Computations  have  been  carried  out  using  the  one 
dimensional  experimental  value  of  unconstrained  angular  change 
along  the  welded  edge  and  its  equivalent  constrained  welding 
moment  as  an  input  to  the  computer  program. 

These  computed  results  are  compared  to  the  particular 

experimental  values.   >. -^  results  indicate  that  an  analytical 

- 

approach  using  fij  *  .^'lement  method  can  predict  the  distor- 
tion phenomena  of  a  complicated  two  dimensional  structure  with 
a  reasonable  degree  of  accuracy. 


Several  recommendations  are  made  concerning  further 
investigations  aimed  at  evaluating  relevant  value  of  angular 
change  and  equivalent  constrained  welding  moment. 
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I   INTRODUCTION 
A.   Background 

Welding  is  used  extensively  in  the  fabrication  of  many 
structures,  including  ships,  airplanes,  buildings,  pressure 
vessels,  etc.,  providing  many  advantages  over  the  techniques 
such  as  riveting,  casting,  and  forging.   However,  distortion 
problems  are  always  encountered  during  welding.   Distortion 
often  prevents  the  ■aCiiiHvenier+  of  design  tolerances,  reduces 
joint  strength  by  mismatching,  and  imparts  initial  deflection 
into  structural  members. 

A  typical  structural  member  is  a  panel  structure,  which 
is  composed  of  a  flat  plate  and  longitudinal  and  transverse 
stiff eners  fillet  welded  to  the  plate,  as  shown  in  Figure  1. 
Angular  changes  produced  along  the  fillet  welds  cause  out-of- 
plane  distortion  of  the  panel.   The  excessive  out-of -plane 

distortion  reduces  the  buckling  strength  of  the  panel  that 

(3) 
is  subjected  to  compressive  loading.      Corrugation  failures 

of  bottom  plating  in  some  welded  cargo  ships  are  believed  to 

be  due  primarily  to  the  reduction  of  buckling  strength  of  the 

plating  with  excessive  initial  distortion.   -  '   ' 

Welding  distortion  is  very  complicated  phenomena 

associated  with  many  far  ors:. 

1.  Heat  conduc'  «  by  of  weldment  which  controls  the 
distribution  of  temperature  gradient  of  welded  metal 

2.  Thermal  expansion  coefficient  which  controls  the 
expansion  and  shrinkage  of  welded  metal. 
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3.  Moduli  of  elasticity  which  governs  the  rigidity 
against  deformation. 

4.  Yielding  strength  which  governs  the  size  of  the 
plastic  zone  near  the  welds. 

5.  Degree  of  constraint  which  governs  the  freedom  of 
movement  during  weld  deformation. 

6.  Amount  of .  heat  input  used  which  governs  the  size  of 
the  melted  zone. 

7.  Amount  of  weldment  used. 

Therefore,  it  is  extremely  difficult  to  analyze  the 
welding  deformation  phenomena  by  a  pure  analytical  way.   Yet, 
only  one  dimensional  deformation  of  panel  structures  has  been 
analyzed  with  the  assumptions  of  elastic  deformation  during 
welding  process.     Obviously,  welding  deformations  are 
neither  pure  elastic  nor  pure  plastic  phenomena,  but  the 
plastic  zone  near  the  weldments  is  so  small  compared  to  panel 
structural  size  that  it  is  believed  to  be  closer  to  the 
elastic  behavior. 

Therefore,  the  analytical  approach  of  the  one  dimensional 
deformation  prediction  can  be  treated  as  an  elsstic  deforma- 
tion by  considering  moduli  of  elasticity,  degree  of  constraint, 
and  amount  of  weldment  used. 

B.   Previous  Investigations 

(1  3) 

Masubuchi,  et.  al.   '    have  investigated  the  one 

dimensional  plate  deformation  during  welding  with  a  minimum 
strain  energy  concept  using  elastic  bending  theory. 
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a..      Free   Joint 


JZlCrlD      BEAM 
/////////////////////////////////////////////////////// 


Bottom  plat£ 


b.   Constrained  Joint 
(framed  structure) 


Figure  2   Distortion Caused  by  Angular  Change  in  Two  Types 
of  Fille-  Welded  Structures 
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If  a  fillet  joint  is  free  from  external  constraint,  the 
joint  simply  bends  to  a  polynomical  form  having  a  knuckle 
at  the  weld  as  shown  in  Figure  2a.   However,  if  the  joint  is 
constrained  by  some  means,  a  different  type  of  distortion  is 
produced.   For  example,  if  the  stiff eners  are  welded  to  a 
rigid  beam,  as  shown  in  Figure  2b,  the  angular  changes  at 
fillet  welds  cause  wavy  or  arc-form  distortion. 

In  the  simplest  case  in  which  sizes  of  all  welds  are  the 
same,  the  distortions  of  all  spans  are  equal  and  the  distor- 
tion, w,  which  can  be  expressed  as  follows: 

$/£  =  [1/4  -  (x/£  -  1/2) 2]  •  9  (1) 

where  6  =  deflection  of  plate 

9   =  angular  change  at  a  fillet  weld  with  constraint 
I   =  length  of  span 

The  amount  of  angular  change,  0.,  in  constrained 

structures  is  smaller  than  that  in  a  free  joint,  6  .   This 

J       o 

indicates  that  a  certain  amount  of  energy  is  necessary  to 

decrease  the  angular  change  from  6   to  8 .   If  the  necessary 

energy  is  represented  by  U  ,  it  may  be  given  by  the  following 

w 

equation. 

e0-e   du 
uw=  <      dTF^eT  d(eo-9)  (2) 

0      o 

On  the  other  hand,  the  strain  energy  stored  in  the 
constrained  plate,  U  ,  can  be  expressed  using  the  elastic 
beam  theory. 

up  =  ^e2  (3) 
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2 


where  E1  =  E/l  -  v 
I  =  t3/12 


Since  U   increases  but  on  the  contrary  U   decreases  as 

p  J      w 

the  constrained  angle  8  increases,  the  condition  for  equili- 
brium of  this  system  requires  that  the  total  strain  energy 
U.  =  U  +  U   should  be  minimum.   Furthermore,  for  the  sim- 
plicity of  the  problem,  the  ratio  of  incremental  welding 
energy  change  to  angular  change  is  assumed  to  be  linear  as  in 
the  following  equation. 

dU 

=  C(0  -9)  (4) 


d(0  -9)     v  o 

where  C  =  a  coefficient  determined  by  the  weight  of  the 
deposited  metal  and  by  the  welding  procedures. 
From  equations  (2)  and  (4) ,  welding  energy  per  unit 
length  of  width  can  be  expressed  in  terms  of  C,  9,  and  9  , 

uw  "  I  (e o-0)2  (5) 

w     2    o 

Accordingly,  the  condition  of  equilibrium  is  as  follows: 

^«  =  C(6o-9)  +  2^1=  0  (6) 

From  equation  (6) ,  relation  between  9  and  9   can  be 


expressed, 

9 
O 


1  +  2E^I/C 


(7) 


The  values  of  C  can  be  determined  by  using  relationship 
obtained  by  experimentally  measuring  9   and  9 .   Empirical 
formula  for  the  calculation  of  C  has  been  proposed  by 
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Figure  3   Variation  of  Angular  Change  of  Free  Fillet  Welds, 
60,  As  a  Function  of  Plate  Thickness,  t,  and 
Weight  of  Electrode  Consumed  per  Weld  Length,  W. 


18 


.  (3) 

Masubuchi. 


C  =  t4/(l+W/5)  kg-mm/mm  radian         (8) 

where  W  =  weight  of  electrode  consumed,  gram/cm. 
t  =  thickness  of  plate,  mm. 

The  experimental  results  of  determining  the  unconstrained 
angular  change  as  a  function  of  plate  thickness,  t (mm) ,  and 

weight  of  electrode  consumed  per  weld  length,  W (gr/cm) ,  have 

(2) 
been  given  by  Hirai ,  et.  al.     as  in  Figure  3. 

Therefore,  for  a  given  plate  geometry  and  welding  con 

condition,  the  deflection  of  the  one  dimensional  case  can 

easily  be  estimated  by  using     equations  (1) ,  (7) ,  (8) , 

and  Figure  3. 

C.   Aim  and  Purpose  of  Present  Studies 

As  stated  previously,  welding  deformation  has  been 
analyzed  only  for  the  one  dimensional  case,  but  in  reality 
most  of  the  structural  members  are  encountered  in  two 
dimensional  deformation,  in  which  the  deformation,  w,  varies 
along  the  x-  and  y-directions  as  well  as  the  unconstrained 
angular  change,  6,  and  equivalent  constrained  welding  moment, 
C.   With  the  idea  of  previous  investigations  which  state  that 
the  total  strain  energy  in  the  plate  during  welding  deforma- 
tion becomes  minimum,  it  is  possible  to  extend  the  one 
dimensional  deformation  analysis  to  the  two  dimensional  case. 
However,  the  simple  mathematical  solution  using  plate  bending 
equation  is  rather  difficult  due  to  the  uncertainty  of  the 
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equivalent  load  and  the  boundary  conditions  relevant  to  the 
welding  deformation  problem. 

The  difficulties  can  be  easily  overcome  by  formulating 
the  problem  into  a  finite  element  model  because  it  can  be 
self -ad justed  to  take  care  of  the  equivalent  load  and  equiva- 
lent sitffness  matrix  due  to  welding.   Furthermore,  in  this 
way  the  relation  of  deformation  phenomena  to  the    constrained 
angular  change,  6,  and  the  equivalent  constrained  welding 
moment  C  can  be  easily  seen. 

Therefore,  it  is  intended  to  handle  the  two  dimensional 
deformation  with  the  finite  element  displacement  technique. 
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II  FORMULATION  OF  FINITE  ELEMENT  EQUATION 

A.   One  Dimensional  Case 

In  this  study,  computer  programs  are  programmed  only 
for  the  two  dimensional  finite  element  analysis,  but  the 
one  dimensional  finite  element  formulations  are  presented  to 
visualize  the  inside  of  the  complex  matrix  assembling  and 
computational  procedure;  furthermore,  to  compare  the  results 
to  that  of  the  previous  investigators. 

For  the  one  dimensional  case,  as  shown  in  Figure  2b,  if 
deflections  are  elastic  and  small,  the  deflection  can  be  cal- 
culated using  elastic  theory  when  the  welding  equivalent  loads 
and  the  boundary  conditions  are  prescribed.   Even  though  the 
deformations  are  assumed  to  be  elastic  for  the  welded  struc- 
tural case,  the  loading  and  boundary  conditions  to  be  applied 
for  this  problem  are  not  clear. 

As  introduced  by  Masubuchi,  et.  al.    ,  if  the  concept 
of  equivalent  constrained  welding  moment,  C(lb-in/in  rad.) 
is  being  used,  the  uncertainty  of  loading  conditions  may  be 
overcome;  for  the  boundary  condition  it  may  be  safe  to  assume 
to  be  simply  supported  for  which  only  the  angular  change  along 
the  welded  edge  can  be  allowed.   Furthermore,  for  simplicity 
of  the  problem,  the  deflections  are  considered  only  for  the 
case  of  one  span  length  of  a  simply  supported  beam  to  elimi- 
nate the  complexity  of  interaction  due  to  the  statically 
indeterminate  structural  effects,  as  in  Figure  4a. 
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"a.   Simplified  one-span  beam  element 


b.   Local  coordinate  and  nodal  points  for  a 
finite  beam  element 


Figure  4   One-span  Beam  Element  and  the  Finite 
Element  Local  Coordinate 
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For  the  one  dimensional  plane  stress,  the  strain  energy 
can  be  expressed  simply  by: 

Up  =  1/2  ///  cxaxdxdydz  (9) 

Using  the  stress-strain  and  strain-curvature  relation  for  the 

beam  element, 

a     =  Ee  (10a) 

x     x 

ex   =  "z  — 2  (10b) 

dx 

Therefore,  the  strain  energy  can  be  expressed  in  terms  of 

deflection  and  rigidity  of  the  plate: 

a 


P         2         Q 


A2 

d  w 


A     2 

dx 


dx  (11) 


where 


t/2      b  b   3 

.1   =   2/        /    z^dzdy  =  ~- 
0        0  ±Z 

On  the  other  hand,  the  welding  energy  can  be  expressed 

in  terms  of  9  ,  0,  and  C,  as  in  equation  (5) .   For  this 

problem  it  is  more  convinient  to  express  the  welding  energy 

in  terms  of  deflection  rather  than  the  angular  changes.   From 

equation  (5) , 

a   b         9 

u,,  =   f   /   |(eo-e)    dydx 

w    0   0 
the  constrained  angular  change  can  be  expressed  by  the 
equation 

9  =  8w/9x  (12) 

Therefore,  the  welding  energy  expressed, 
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U 


w 


be 
2 


a 

/  0 
0 


3w 

dx 


)      dx 


(13) 


where  the  values  of  C  and  0   are  constant  which  will  be 

o 

determined  by  the  particular  welding  condition  and  the 

geometry. 

The  equilibrium  condition  for  the  system  requires  that 

the  total  strain  energy  has  to  be  minimized f  in  other  words, 

the  variation  of  the  total  strain  energy  has  to  be  zero, 

which  becomes 

<$Um  =  6U   +  SU   =0  (14) 

T     p     w 

using  equations  (11)  and  (13) ,  equation  (14)  becomes 

2 


<SUT  =  6 


EI 
2 


a 
/ 
0 


A2 

d  w 


dx' 


dx 


+  6 


be  r    ,Q    dw,2 
o—  /  (9  -  -g— )   dx 
2   o    °   dx 


But  the  variation  of  the  constant  term  is  zero,  therefore, 

„T  f    ,d  6wx  ,d  w,  ,   ,  be  ,  ,d6wN  ,dwx  ,     ,   r    ,d6wN  ,   ,nc-v 

EI  {  ("^^)(dx^)dx  +  T  [  {~^]  (dx")dx  -  bc  f0  (-dir)dx  (15) 

For  a  beam  element  as  shown  in  Figure  4b,  the  deflection 

can  be  assumed  to  be  a  general  cubic  function  of  x, 

2       3 
w  =  c,  +  c„x  +  c^x   +  c.x  (16) 

and  when  the  nodal  displacements  and  the  nodal  angles  are 
the  unknown  variables  which  are  to  be  calcualted,  the  coeffi- 
cients (c, ,  c2,  c_,  c.)  in  equation  (16)  can  exclusively  be 
expressed  in  terms  of  the  nodal  displacements  and  angles 
considering  the  requirements  of  continuity  at  the  nodal  points 


w(x=0)  =  w,  =  c. 
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^(x=0)  =  wXi  -  c2 

2       3 
w(x=a)  =  w?  =  c,  +  c~a  +  c_a   +  c.a 

3w ,    »  ~       -.2 

■-—  (x=a)  =  V7   =  c0  +  2c..a  +  3c.a 

oX  X~      z       J        4 

The  above  requirements  can  be  expressed  in  the  matrix  form, 
q  =  T  C  (17) 

where 

"10   0    0 


0  10    0 

1  a   a    a 


0   1   2a   3a' 


(17a) 


e         \ 

ci 


c  = 


f*. 


(17b)      and 


W- 


W 


x. 


w. 


w 


(17c) 


From  equation  (17) ,  C  can  be  expressed 
C  =  T_1  q 


(18) 


If  the  new  matrix  G  is  defined  such  that 
G  =  T"1 


(19) 


then, 


C  =  G  q 


(20) 
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where 


G  = 


1 

0 

0 

0 

0 

1 

0 

0 

■>/    2 
3/a 

-2/a 

->,    2 
3/a 

-1/a 

2/aJ    1/a    -2/a" 


1/a' 


(21) 


The  second  derivative  of  deflection  function  becomes, 

A2 

d  w   0    .  r 
— y  =  2c,.  +  6c.  x 

dx^      J      4 


which  is  in  matrix  form 


A2 

d  w 

a    2 
dx 


=  H  C 


where  H  is  the  matrix  defined  by 
H  =  (0   0   2   6x) 

then, 

r„2r  ^  rj2  ^ 

6(H  C)  (H  C)  = 


d    Sw" 

/i    2 
dx 

\                      J 

fA2      \ 

d  w 

A     2 

dx 

(6CT) (HTH  C) 


(22) 


(23) 


but,  using  the  equation  (20)  , 


d    Sw 

A     2 

(dx      J 

d  w 

A     2 

[dx    J 

=  (6qT) (GT  HT  H  G  q) 


(24) 


The  matrix  expression  for  the  first  derivative  of  the 
deflection  w  can  be  expressed  in  the  same  manner, 

2 


dw 
dx 


c2  +  2c-.x  +  3c.x 


using  F_  and  C. 
dw 


dx 


=  F  C 


(25) 
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where 

P  =  (0   1   2x   3x2) 
Finally,  using  equation  (2  0) 

and 


m      mm 

=  (6q  )  (GT  Fx  F  G  q) 


(26) 


(27) 


d6w' 
[dx    J 

'dw 
,dxl 

(28) 


From  equations  (24) ,  (27)  and  (28) ,  the  minimum  energy 
equation  (15)  can  be  written  in  the  matrix  form, 


EI  /  (<5qT)  (GT  HT  K  G  q)dx 
0    ~      ~ 


+bc  /  (<5qT) (GT  FT  F  G  q) dx 


-be  /  (<5qT)  (GT  FT)dx  =  0 
0 


(29) 


T 
Noticing  that  the  matrix  q,  q  ,  G,  and  G,   is  not  the 

function  of  x,  but  only  of  the  function  of  dimension  of  the 
plate,  and  displacement  and  angle  at  the  nodal  points;  there- 
fore, equation  (29)  can  be  simplified  to 


6qT  (EI  GT  [/  (HT  H)dx]G  q 


+bc  GT  [ /  (FT  F)dx]G  q 

~   Q 

-be  GT  [/  (FT)dx]}  =  0 

0  "* 


(30) 
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in  which  the  only  integrations  necessary  to  be  carried  out 
are  the  bracketed  quantities,  which  become  after  integration, 


/(H1  H)dx  = 
0  ~  ~ 


0  0  0 
0  0  0 
0   0   4a 


0 
0 
6a' 


0   0   6a2   12a3 


(31) 


/(Fx  F)dx 

0  ~     ~ 


0   1 


2a    3a" 


0   2a    4a2   6a3 

n   t  2  a    3   n  4 
0   3a    6a    9a 


(32) 


/(Fx)dx  = 

0  ~ 


0 
1 

2a 
3a: 


(33) 


Equation  (30)  is  the  final  equation  for  one  element 
formulated  by  the  finite  element  displacement  method.   The 
first  quantity  is  defined  as  a  stiffness  matrix  for  an 
element,      the  second  is  the  additional  stiffness  matrix  due 
to  welding,  and  the  third  is  the  equivalent  load  due  to  weld- 
ing.  They  are  as  follows: 
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Stiffness  matrix  for  an  element, 

k  =  EI  GT[/(ET  K)dx]G 
~  0 


=  EI 


1             3 
12/aJ 

c/    2 

6/a 

-12/a3 

en    2 

6/a 

6/a2 

4/a 

cj    2 
-6/a 

2/a 

-12/a3 

o    2 
-6/a 

12/a3 

-6/a2 

>   6/a2 

2/a 

a/    2 

-6/a 

4/a  j 

(34) 


Additional  stiffness  matrix  due  to  welding, 
kTT  =  be  GT[/(FT  F)dx]G 


0  0  0  0 

0  0  0  0 

0  0  0  0 

0  0  0  be 


(35) 


Load  matrix  due  to  welding, 

T  a   T 

Q  =  be  G1 [/  F1  dx] 


0 

0 

bc6 


(36) 


At  this  stage,  assembling  of  equation  (3  0)  has  to  be 
considered  when  the  one-span  beam  is  divided  by  many  finite 
elements,  but  the  one-span  beam  element  of  Figure  4a  is 
chosen  to  be  two  finite  elements  as  in  Figure  4b;  then,  only 
element  1  can  be  used  to  predict  the  deflection  and  angular 
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change  at  the  nodal  points  by  symmetric  argument.   In  this 
case,  the  simultaneous  equations  to  be  solved  are  (from 
equation  (30) ) , 

6wl[^T  (l  wl  +  wxn  "  I  w2  +  Wx>J  "  ° 
a  1  2 

'VtF    CIw1   +    2%   "   lw2   +wx2)]    =    ° 

«„2  [M|  (_|       .         +  a       _  „  0 

a  12 

*V^  (|  wx  +  w^  -  §  w2  +   (2  +  ||£)wX2  -  bc6o)]   =  0 

(37a-d) 

From  the  simply  supported  boundary  condition  at  node  2, 

w~  =* 0;  therefore,  6w~  =  0. 

From  symmetry  condition  at  node  1,  w   =0;  therefore, 

xl 
6w   =0. 
xl 

Furthermore,  the  displacement  at  node  1  and  the  angle  at 
node  2  is  not  known,  which  means  that  the  variations  of  these 
variables  are  not  zero;  therefore,  the  only  way  to  satisfy 
equation  (37a)  and  (37d)  is  the  bracketed  value  in  each  equa- 
tion has  to  be  zero,  which  leads  to  final  equation  for  w 

and  w   , 
X2 

2 

—  w,  +  w   =0 
a   1    x2 

I  wi +  (2  +  §!>wx2  ■  °  (38a"b) 
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From  equation  (3  8a-b) ,  displacement  at  the  middle  point  and 

angular  change  at  the  edge  expressed  as, 

wl  "  2(1  +  El/abc)  ^*> 

w    =  *m — .  r^n-  i    w  ^  (40) 

x     2  (1  +  El/abc) 

It  is  worth  to  compare  the  results  of  equations  (6)  and 
(7)  with  equations  (39)  and  (40) .   In  equations  (39)  and  (40) , 

EI   =  bt3/12  :  t3/12 
abc     abc       ac 

and  a  =  1/2   is  to  be  used.   Then, 

Wl  =  4(1  +  2EI/JLC)  (   ' 

•wx2   "   1   T   2ElAc  <42> 

where  I  is  defined  by  unit  width  basis,   '    which  is 
I  =  t3/12. 

Equation  (42)  agrees  with  equation  (7) ,  as  also  does 
equation  (41),  if  x  =  1/2   and  6   of  equation  (7)  is  plugged 
into  equation  (1) .   The  negative  sign  of  equation  (41)  means 
that  the  coordinate  and  the  positive  angle  at  the  edge  has  to 
be  chosen,  as  in  Figure  4a. 

As  shown,  the  results  from  the  finite  element  method  is 

(1  3) 
exactly  the  same  expression  proposed  by  Masubuchi ,  et.al.   ' 

for  the  one  dimensional  deformation,  and  in  this  way,  the 

equivalent  loads  and  the  stiffness  matrix  due  to  welding 

become  easily  visualized. 
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B.   Two  Dimensional  Case 

1.   Formulation  of  the  Problem 

In  the  two  dimensional  analysis,  all  procedures  are  the 
same  as  that  of  the  one  dimensional  analysis  except  that  the 
integrations  are  carried  out  by  the  numerical  methods  and  the 
methods  of  assembling  are  mentioned.   The  deflection  may  be 
approximated  by  the  bi-cubic  function  for  one  element,    ' 

as  in  Figure  5. 

2       2       3 

w  =  c,  +  CpX  +  c,y  +  c.xy  +  c^x  +  cfiy  +  c-x 

2        2        3       3        2  2 
+  c8x  y  +  cgxy   +  c1()y   +  c-^x  y  +  c12x  y 

+  ci3xy    +  ci4x  y    +  ci5x  y    +  cigx  y  (43a) 


9~x  =  C2  +  C4Y  +  2c5X  +  3c7X   +  2c8xy  +  C9Y 

+  3c11x  y  +  2c12xy   +  c13y   +  3c14x  y 

3         2  3 
+  2c, rxy   +  3c, ^x  y  (43b) 

lo         lb 


|~  =  c3  +  c4x  +  2c6y  +  c8x2  +  2c9xy  +  3c1Qy2 

3         2  2         3 

+  c, ,  x   +  2c,  2X  Y  +  3c,  ~xy   +  2c,  .x  y 

2  2         3  2 
'  +  3c.  rx  y   +  3c,  fix  y  (43c) 

32w  2 

=  c.  +  2cpx  +  2cQy  +  3cnix   +  4c, 0xy 


8x8y    "4    "^8"        ^9*  '  """irk    ""12 

2         2  2 

:13y  +  6c14x  y  +  6c15xy   +  9c16 


2         2  2         2  2 

+  3c,  ~y  +  6c,  .x  y  +  6c,  ^xy   +  9c,  ^.x  y        (43d) 
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0  15  -  1^"h4- 


4 


4—*% 


(o,o; 

^4=  tJxMI 


(a,o) 


Figure  5   Local  Coordinate  of  a  Two  dimensional  Rectangular 

Finite  Element  and  Generalized  Displacement  Numbering 
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From  equation  (4  3a-d) ,  the  matrix  equation  can  be  defined 

such  that  at  each  nodal  point,  the  generalized  displacements 

(deflections,  angles  in  x  and  y,  and  twisting  angles)  are 

expressed  in  terms  of  the  geometry  of  the  element  and  the  some 

unknown  coefficients  C,  which  is  the  form  of, 

q  =  T  C  (44) 

where  T  is  16  by  16  matrix,  C  is  16  by  1  matrix,  and  q  is 
16  by  1  matrix  defined  as  in  Appendix  A. 

From  equation  (44),  the  unknown  coefficient  matrix,  C, 
can  be  expressed  in  terms  of  geometry  and  the  generalized 
displacements , 

C  =  Gg  (45) 

where 

.  G  =  T_1  (46) 

However,  the  actual  calculation  of  the  inverse  of  T  is 

carried  out  by  the  computer  in  this  analysis. 

(8  9  10) 
From  the  elastic  plate  theory,   '  '     the  strain  energy 

stored  in  the  plate  during  the  elastic  deformation  is  ex- 
pressed by 


D 


a  b 


p   2  0  0 


3  w    3  w 


-  2(l-v) 


Jx     3y 

This  can  be  written  in  the  matrix  form, 
1   a  b 

U   =  =■  /  /  K   D  K  dydx 
p      ~  *   ~  ~  ~* 

where 


^2   ~2 

3  w.  3  w 

2    2 

3x^  dy 


32w 


3x3y 


dxdy 


(7) 


which  is 


0  0 


(47) 
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K  =  - 


(      2      ) 
3   w 

s    2 
3x 

*2 
3   w 

a    2 
3y 

032w 

■"3x3y 

-H  C 


in  which  H  is  given  by  Appendix  A, 


f. 


D'  = 


Ef 


12  (1-v  ) 


0 

0 

1-v 


=  B 


in  which 

3  - 

and  * 


Ef 


12  (1-v  ) 


D  = 


1 

V 

0 

V 

1 

0 

0 

0 

1-v 

(48) 


(49) 


Therefore,  equation  (4  8)  becomes  after  taking  the 

variation  with  respect  to  generalized  displacement, 

a  b 
6U   =  /  f    (6K  ) 3  D  K  dydx 
P    0  0 

a  b    _,    _ 
-  /  /  (6CT) (HT  6  D  H  C)dydx 
0  0 


T    T  a  b  _ 
6q   B  G  [/  /  H   D  H  dydx]G  q 
0  0  ~ 


(50) 
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The  welding  energy  for  the  elements  shown  in  Figure  6b 
can  be  written  as, 


a  C 


b  C 


uw  =  {  -i  VV  dx  +  {  "I  (W  dy 


Therefore, 


U. 


C   a 

^y 


W 


T 


/  (e~  - 


3wx  ^    ,  cx  ,  ,n 


9_w 

ax 


)  dy 


which  becomes, 


uw  = 


,8v>  2         ,    cx    r     ,9w.  ^ , 
o    /     (^7)     dx    +   T  /     (3^7)     dY 


2 


9y' 


9xJ 


(51) 


C8      /|^dx  +   C0      /   22  dy 

y  o  0  3y  xo08xT 


*r(c     +  c  )e2 

2      y  x      o 


(52) 


and  taking  the  variation, 


w 


a 
C      / 

_y  0 

9Sw 

I  3yJ 

(                \ 

9w 

UyJ 

dx  +  C   / 
X  0 


{    9xJ 

'aw' 
[syj 

dy 

c  e  / 
y  °  n 


9<Sw 


isyj 


dx  +  C 


x  o 


a 


9_6w 
9x 


dy 


(53) 


in  which  the  first  bracketed  term  is  the  additional  stiffness 
matrix  due  to  welding,  and  the  second  is  equivalent  load 
matrix  for  the  welding.   In  matrix  form,  equation  (53)  can  be 
written,  if  some  matrices  defined  as  below: 


^  =  F   c   =  cT  FT 
9x   £y  ~    —  ~y 


(54a) 
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3w  T   T 

22L  =  F   c  =  C   F 

3y   ~x  ~   ~     — x 


(54b) 


where  the  subscripts  x  and  y  denote  that  the  integration 
should  be  carried  out  along  the  x-  and  y-direction,  respec- 
tively, and  the  matrix  F   and  F   are  defined  in  Appendix  A. 
Therefore, 


'3<5w^ 

=  (6CT)  (FT  F   C) 
*—  '  v«~y  ~y  ~/ 


(6qT)GT  FT  F   G  q     (55a) 
x.  —  ~y  -~y  '—  -3> 


3  6w 
1     3xi 

3w 
,3x 

k.             J 

(6CT) (FT  F   C)  =  (6qT)GT  FT  F   G  q     (55b) 

~    —X  —X  ~         ~  ~-      ~. X  ~X  —  ~ 


Using  equations  (54a-b)  and  (55a-b) ,  the  welding  term 
can  be  written  as: 


5Uw=  6%     Cy  <£    V     (£x^x)dx]S.l 


+    6qT   C      GT    [/     (FT   F    )dy]G   q 
x    ~~       «     >~y  >--y  ^-  i. 


0 


-    <5qT   C      6      GT[/     (FT)dx] 
~       y     o   -~     rt     —x 


-    6qT   C      9o    GT[/     (F^)dyJ 


(56) 


As  in  the  one  dimensional  case  already  shown,  the  total 

variation  should  be  zero  to  be  in  equilibrium  of  this  system, 

Therefore, 

6Um  =  6U   +  SU   =0 
T      p     w 

which  is  from  equations  (50)  and  (56) , 
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a   b 


tTtt  aT 

Ug   )    G    [3    /    /   H"    DH   dydx   +   C      /      F      F      dx 

~     r         0  0  ~    ~  y  Q    ~x  ~x 

b  a 

+  C   /  FT  F   dy]G  q  -  GT[C  0   /  FT  dx 
^  0  ~Y  ~Y  *   ~  ~        ~   Y  o  Q  -x 

+  C  6   /  FT  dyj)  =  0  (57) 

x  o  Q  -y  *J  J 

These  are  the  16  by  16  simultaneous  equations  which  have 
to  be  solved  if  only  one  element  is  concerned  as  shown  in  the 
one  dimensional  case,  with  the  proper  boundary  and  symmetry 
conditions.   From  equation  (57),  the  stiffness  matrix  for  a 
plate  element, 

k   =   GTg//HTDH   dydx  G  (58a) 

~        0  0  ~    ~ 

The  additional  stiffness  matrix  due  to  welding  along  the 

x-direction  of  the  plate, 

k    =  GT  C   /  FT  F   dx  G       •  (58b) 

-wx   —   y  ~  —x  —x    ~ 

The  additional  stiffness  matrix  due  to  welding  along  the 
y-direction  of  the  plate, 

k    =  GT  C   /  FT  F   dy  G  (58c) 

~wy   ~   x  Q  -y  ~y   J  ~ 

The  equivalent  load  matrix  due  to  welding  along  the 
x-direction  of  the  plate, 

m  ^    m 

Q  •  =  GT  C  9   /  FT  dx  (58d) 

~*   ~   y  0  Q  -x 

The  equivalent  load  matrix  due  to  welding  along  the 
y-direction  of  the  plate, 

Q   =  GT  C  9   /  FT  dy  (58e) 

Xy    ~   x  o  0  ~y 
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*X 
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Figure  6  Various  Finite  Element  Having  Welded  Edge 
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2.  Assembling  Procedure 

For  all  the  finite  elements  except  those  having  the 
welded  edge,  the  stiffness  matrix  necessary  to  be  assembled 
is  the  expression  given  by  equation  (58b) ,  and  for  the  finite 
elements  having  welded  edges,  three  cases  are  involved.   The 
first  case  is  elements  having  welded  edges  along  the  x-direction 
as  in  Figure  6b;  the  second  case  is  along  the  y-direction,  as 
in  Figure  6c;  and  the  tluLrd  case  is  along  x-  and  y-directions 
as  in  Figure  6b.   The  assertibTim)  procedure  is  the  expression 

of  equation  (57)  for  the  overall  plate  which  is  to  be  solved, 

(13) 
which  is  described  in  FEABL  user's  manual.       Also,  from 

the  symmetry  condition,  which  will  be  described  in  the  follow- 
ing section,  only  one-quarter  part  of  the  plate  is  actually 
being  considered. 

Example  of  the  assembling  procedures  for  the  two  finite 
elements  is  given  in  Appendix  B. 

3.  Symmetry  and  Boundary  Conditions 

After  assembling,  relevant  symmetry  and  boundary 
conditions  have  to  be  applied  to  solve  the  simultaneous 
equations,  as  shown  in  the  one  dimensional  case.   Along  the 
welded  edge  of  x  and  y,  the  simply  supported  edge  condition 
requires  that  all  the  displacements  at  that  nodal  point  have 
to  be  zero,  and  for  the  x-directed  edge,  the  angular  changes 
in  the  x~direction  are  not  allowed,   In  just  the  same  way, 
for  y-directed  edge,  the  angular  changes  in  y  are  not  allowed. 
The  symmetry  condition  requires  that  the  first  derivative 
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along  the  symmetric  axis  should  be  zero.   For  example,  the 

symmetry  and  boundary  conditions  for  the  case  given  by 

Figure  7  are  as  follows: 

Boundary  conditions : 

At  node  6,  12,  18,  24,  30,  36,  42,  48,  54 

w  =  0 

w  =  0 
x 

At  node  49-   SQ$    51,  52,  53,  54 
w  =  0 
wy  =  0 

Symmetry  conditions: 

At  node  1 ,  2 ,  3 ,  4 ,  5 ,  6 
w  =  0 

y 

At  node  1,  7,  13,  19,  25,  31,  37,  43 

Wx  =  ° 

At  node  1,  6,  4  9 

w   =0. 
xy 

4 .   Numerical  Integration  Method 

Usually  in  the  finite  element  analysis,  the  integrations 
involved  are  carried  out  by  the  numerical  methods ,  which  can 
be  easily  calculated  using  the  computer.   For  the  present 
analysis,  the  nine  Gaussian  points  are  chosen  for  the  rectan- 
gular finite  element,         as  shown  in  Figure  8c.   In  actual 
integration,  it  is  very  convenient  to  use  the  normalized  coor- 
dinate system,  as  will  be  described  below.   From  Figure  8a  and 
8b,  normalized  coordinate  system  can  be  defined  as 


x  =  |(X  +  1) 


Y  =  fe(Y  +  1) 
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(59a) 

(59b) 

(60a) 

(60b) 

and 

dx  =  j   dX 

dy  =  |  dY 

Therefore,  the  arbitrary  integration 

a  b 
f(a,b)  =  /  /  P(x,y)dxdy  (61) 

0  0 

becomes  in  the  normalized  coordinates, 

,11  , 

f(a,b)  =  ~      f      /  P[|(X+l)f  §(Y+1)]  (62) 

*   -1  -1    A  * 

< Equation  (62)  can  be  integrated  by  numerically  using 
9-Gaussian  points,  which  become 

,33 

f(a,b)  =  2%     I    Z   W.W.  P(x,y)  (63) 

4  i=l  j=l   x  -1 


where 


x  =  |(XjL  +  1)  (63a) 

y  =  |(Y.  +  1)  (63b) 

Wi:      Wl   =   5//9,  W2   =    8//9/    W3   =    5//9  (63c) 

W.:      Wx   =    5/9,  W2   =    8/9,    W3   =    5/9  (63d) 


X, 


i:      X±   =   -/15/5,    X2   =    0,    X3   =    /15/5  (63e) 


:.:      Yx    =    -/I5/5,    Y2   =    0,    Y3   =    /15/5 


(63f) 
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Figure  8   Ordinally  and  Normalized  Coordinates  for  an 
Element,  and  9  Gausian  Points  for 
Numerical  Integration 
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Therefore,  using  the  concept  of  normalizing  and 
numerical  integration  described  above,  the  integrations 
involved  in  the  present  analysis  such  as  equation  (58a-e) 
can  easily  be  written,  which  are  programmed  in  Appendix  D. 

From  equation  (58a) 

3    3 
k  -  ffl  GT  [  Z    Z   W.W.  P(x,y)]G  (64) 

~   i=l  j=l  1  : 


where 


P(x,y)  =  HT  D  H 
x,  y,  W. ,  W.,  X.  and  Y.  are  defined  by  equation  (63a-f ) . 

From  equation  (58b) , 

aC       3 
kwv  ■  -o2-  G  I  Z   Wi  Pv(x,y)]G  (65) 

~WX      Z     ■ —      n    1  ^~X        ^ 

1=1 

where 

P  (x,y)  =  FT  F 
~x   ■*    ~x  ~x 

x  =  |(Xi  +  1) 

y  =  b 

W.  and  X,  are  defined  by  equation  (63c)  and  (63e) ,  respectively 

From  equation  (58c) , 

bC  3 

*wy   =  -1T£    [4  Wj    Py(x'Y)]^  (66) 


where 


P    (x,y)    =   FT   F 
— y        ' J  '  ~y  ~y 


X    =    a 

Y  =  ICYj    +   1) 


45 

W.  and  X.  are  defined  by  equations  (63d)  and  (63f ) , 
respectively. 

From  equation  (58d) , 

r=l 
where 

P0  Cx,y)  =  FT 

x  =  |(X±  +  1) 

y  =  b 

W^  and  X^  are  defined  by  equation  (63c)  and  (63  ) , 

respectively. 

From  equation  (58e) , 

b0  C       3 
Qy  =  — 1^-  GX[    E  W.  P£y(x,y)]G  (68) 

where 

x  =  a 

y  =  |(Y.  +  1) 

W.  and  Y.  are  defined  by  equations  (63d)  and  (63f ) . 

As  noticed  above,  numerical  integration  using  the 
normalized  coordinate  can  simplified  the  computing  procedure 
by  factoring  out  a  and  b. 
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5.   Computer  Programs 

Computations  of  the  ordinally,  additional  welding,  and 
equivalent  welding  load  matrix  of  equation  (58a-e)  have  been 
programmed  by  using  the  numerical  integration  technique  with 

normalized  coordinates  as  listed  in  Appendix  D.   In  addition 

(13) 
to  these,  FEABL  Programs     in  assembling  of  the  relevant 

matrices  and  solving  the  simultaneous  equations  with  proper 

boundary  and  symmetry  conditions  are  also  used  and  listed  in 

Appendix  D. 

The  Programs  listed  in  Appendix  D  are  only  relevant  to 
the  panel  structure  welding  with  rectangular  finite  elements. 
Further,  the  nodal  numbering  system  should  follow  the  way 
given  by  the  sample  analysis,  which  means  that  the  numbering 
should  start  from  origin  to  y-direction,  as  shown  in  Figure  7. 

Input  data  for  this  analysis  are  the  overall  dimensions 
of  one  panel  structure  in  the  x-  and  y-direction,  thickness 
of  the  plate,  Young's  modulus  of  elasticity,  Poisson's  ratio, 
number  of  elements  divided  in  x~  and  y-direction,  equivalent 
welding  moment  in  x-  and  y-direction,  and  free  joint  angular 
change  at  the  edge.   The  input  data  cards  are  described  in 
Appendix  C. 

Outputs  are  the  displacements,  angular  changes  in  x-  and 
y-directions ,  and  the  twisting  angles  at  each  nodal  point,  as 
shown  in  Appendix  E. 

One  thing  to  be  noticed  here  is  that  for  a  quarter  part 
of  the  plate,  as  shown  in  Figure  9,  the  proper  boundary  and 
symmetry  conditions  at  the  nodal  points  are  to  be  generated 
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by  subroutine  HOLD  given  in  Appendix  D;  therefore,  as  long  as 
boundary  conditions  of  simply  supported  and  symmetry  conditions 
as  described  are  concerned,  it  is  not  necessary  to  be  consi- 
dered in  the  input.   Therefore,  no  matter  how  the  number  of 
finite  elements  are  to  be  chosen,  the  only  input  data  will  be 
as  described  in  equation  (69a-h) .   Furthermore,  to  use  this 
program,  the  coordinate  and  a  quarter  part  of  the  plate  should 
be  chosen  as  the  same  way  as  shown  in  Figure  9. 

The  units  used  in  the  sample  program  are  pounds,  inches, 
and  radians.   As  results,  for  the  output  data,  the  deflections 
are  in  inches,  the  first  derivatives  (angular  changes)  are  in 
radians,  and  the  second  derivatives  with  respect  to  x  and  y 
(twisting  angles)  are  also  in  radians. 

•  Program  language  used  is  FORTRAN  IV,  and  the  computations 
are  carried  out  by  an  IBM  370/155  system. 
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III   RESULTS 

As  computer  inputs,  one  steel  panel  structural  dimension, 

as  shown  in  Figure  10,  for  which  the  experiment  conducted  and 

(4  19) 
weight  of  consumed  welding  rods  are  given,  3  gr/cm. , 

are  used  to  predict  the  deformations,  angular  changes,  and 

twisting  angles  at  the  nodal  points.   Here,  the  value  of  the 

one  dimensional  free  joint  angular  change  along  the  welded 

edge,  0  ,  has  been  read  fzr~  the  results  of  Figure  3  and  the 

equivalent  constrained  welding  moment,  C,  is  being  evaluated 

using  equation  (8) .   The  numerical  values  of  the  input  data 

are  as  follows: 

Half-length  of  the  plate  in  the  x-  and  y-directions , 
respectively,  L  =  16  inches  and  L   =10  inches. 

Number  of  finite  elements  in  the  x-  and  y-directions, 
respectively,  M  =  8  and  N  =  5. 

Equivalent  constrained  welding  moment  from  equation  (8) 

r  t4     _  (25.4  x  3/8)4  _     .  kg-mm 

1  +  W/5       1+3/5     '  DZUU  mm  rad. 

where  W  is  given  by  3  gr/cm  from  experiment. 
Therefore, 

C  =  1.13  x  10   lb-in/in  rad. 
Poisson's  ratio: 

v  =  0.3 
Young's  modulus: 

E  =  30  x  io6  psi. 
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Thickness  of  the  plate: 

t  =  3/8  inch.  (69  a-h) 

The  results  of  the  analysis  using  the  above  data  are 
listed  in  Appendix  E,  and  these  results  are  plotted  in  the 
Figures  11  and  12  to  be  compared  with  experimental  results. 
As  can  be  noticed  from  Figures  11  and  12,  the  finite  element 
results  are  lower  value  than  the  experiments,  and  the  maximum 
deflection  at  the  midpoint  of  the  plate  which  is  the  most 
important  deflection  in  reality  is  only  a  half  of  that  from 
the  experiment. 

To  visualize  the  behavior  of  the  deformations  as 

functions  of  6   and  C,  series  of  computations  are  conducted 

with  different  combinations  of  6   and  C  as  shown  in  Table  1. 

o 

Yet  remaining,  the  rest  of  the  input  data  are  unchanged,  and 
the  maximum  deflections  at  nodal  point  1  are  plotted  in  terms 
of  6Q  and  C  as  in  Figures  13  and  14,  and  the  maximum  angular 
changes  at  node  6  and  4  9  are  plotted  in  the  Figures  15  and  16. 
Only  the  maximum  deflections  at  node  1,  and  the  maximum 
angular  changes  at  node  6  and  4  9  from  the  results  of  series 

computations  are  tabulated  in  Tables  2  and  3. 

-3 
Finally,  computed  results  using  9   =  55  *  10    and 

C  =  4.2  x  10   are  listed  in  Appendix  E,  and  also  plotted  in 

Figures  11  and  12  to  compare  with  the  experiment.   As  can  be 

seen  from  Figures  11  and  12,  these  results  are  very  close  to 

the  experiment  and  the  reason  of  choosing  the  above  9   and  C 

will  be  discussed  in  the  following  section. 
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IV   DISCUSSION  OF  RESULTS 

As  has  been  mentioned  in  the  previous  section,  the 

4  -3 

results  using  C  =  1.13  x  10   and  6   =  45  x  10    are  only  a 

half  of  the  experimentally  measured  values  as  shown  in 
Figures  11  and  12.   The  reasons  of  this  difference  are  not 
yet  known.   These  may  result  from  the  assumption  of  the  elas- 
tic behavior  of  the  welded  deformation,  or  from  the  uncer- 
tainty of  the  experimental  measurement  itself  because  the 

deflections  and  angular  changes  which  are  only  of  the  order 

-1  -2 

of  magnitude  of  10    inches  and  10    radians,  respectively. 

Yet,  other  possibilities  are  from  the  value  of  Q      and  C. 
c  o 

To  investigate  the  behavior  of  the  deformation  with 
respect  to  80  and  C,  series  of  computations  have  been  con- 
ducted.  The  plotting  of  maximum  deflections  at  node  1 
versus  C  with  changing  8   as  in  Figure  13  shows  some  insight 
of  deflection  behavior  with  varying  C.   In  here,  it  is  ob- 
vious that  the  deflections  are  more  rapidly  changed  with 
changing  0   than  C  changes,  which  means  that  the  deflections 
are  more  strongly  dependent  function  of  9   rather  than  C. 
Furthermore,  the  deflection  reaches  rapidly  to  an  asymptotic 
value  as  C  increases,  which  can  be  seen  clearly  from  the  one 
dimensional  deflection,  equation  (41)  where  the  term  2EI/ £c 
becomes  small  compared  to  1  as  C  increases;  therefore,  in  the 
limit  case,  the  deflection  can  be  expressed  only  by 

w  =  id    /A. 
o 
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This  behavior  of  the  deflection  with  changing  C  are 

clearly  shown  in  Figure  13.   Also  from  Figure  14,  in  which  the 

maximum  deflections  are  plotted  with  respect  to  6  ,  it  is 

obvious  that  the  deflections  are  rapidly  increased  as  6 

increases. 

If  assumed  that  the  difference  between  experiment  and 

computed  results  are  caused  from  the  possibility  of  different 

values  of  6   and  C,  then  the  question  is  how  to  choose  6   and 

O         '  -a  0 

C,  such  that  the  computed  maximum  deflection  becomes  the  same 
to  that  from  the  experiment.   From  Figures  12  and  13,  the  com- 
binations of  6   and  C  to  satisfy  the  given  value  of  maximum 
deflection  from  the  experiment  can  be  predicted.   However, 
the  best  combination  is  not  known,  but  for  the  known  experi- 
ment value  of  maximum  deflection  it  may  be  safe  to  choose  9 
and  C  such  that  the  changes  of  both  9   and  C  from  the  cal- 
culated value  be  minimized.   In  other  words,  superposing 
Figures  13  and  14  as  shown  in  Figure  17,  the  point  which 

gives  the  same  reading  of  9   and  C  in  both  ways  as  shown  by 

-3 
the  small  circle,  which  is  read  as  9   =  55  x  10    and 

o 
4 
C  =  4.2  x  10  ,  will  be  the  optimal  combination  to  satisfy 

the  given  maximum  deflection.    The  comparison  of  the  results 

using  this  9   and  C  to  the  experiments ,  as  in  Figures  11  and 

12,  shows  that  this  approach  gives  reasonably  good  results. 

As  has  been  seen  from  above,  there  exists  many 

uncertainties  for  the  two  dimensional  welding  analysis,  but 

the  welding  deformation  phenomena  can  be  predicted  with  a 

reasonable  degree  of  accuracy  by  the  finite  element  analysis. 
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Furthermore,  the  empirical  formula  of  C  given  by  equation  (8) 
may  not  be  simply  applied  to  the  two  dimensional  case,  and 
also  does  the  value  of  0  .   If  it  is  so,  measuring  only  the 
maximum  deflection  in  future  experiments,  the  best  combination 

of  6   and  C  can  be  predicted  as  discussed. 

o  c 
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V   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   Conclusions 

From  the  discussion  of  the  previous  section,  it  can  be 
concluded  as : 

1.  The  finite  element  analysis  for  the  two  dimensional 
welding  deformations  which  are  assumed  to  be  elas- 
tic deformation  gives  reasonable  prediction  of 
deflections . 

2.  The  searching  technique  of  optimal  combination  of 

6   and  C  for  a  given  deflection  can  be  used  as  shown 

o  J 

already,  if  there  exists  the  uncertainty  of  6   and  C. 

3.  In  the  future  experiments,  this  analysis  gives  the 
order  of  magnitude  feeling  of  the  deflection  mea- 
sured after  welding.   It  is  extremely  difficult  to 
measure  correct  deflections  in  experiments  because 
the  order  of  magnitude  to  be  measured  is  very  small; 
therefore,  the  inherent  uncertainties  always  exist. 

4 .  The  magnitude  of  importance  of  the  welding  deformation 

variables  such  as  0   and  C  are  clearly  seen,  there- 

o  J 

fore,  more  care  should  be  concentrated  in  measuring 
9   rather  than  C. 

5.  Simple  empirical  formula  for  C  given  by  equation  (8) 
may  not  cover  the  wide  range  of  the  values  of  C; 
furthermore,  as  can  be  derived  from  equation  (41)  it 
may  be  logical  to  express  C  in  terms  of  maximum  de- 
flection, free  joint  angular  change,  and  rigidity, 
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such  as : 

C  =  6  /4w     -  2EI/£  (70) 

o   max 

in  which  6   is  a  function  of  a  given  welding 

condition  and  w    is  the  constrained  maximum  deflec- 
max 

tion  at  midpoint  of  the  span  for  the  same  condition 
of  the  welding  as  of  0  .   Further,  the  elastic 
rigidity  and  plate  dimension  terms  are  expressed  by 
the  term  of  2EI/JI. 


B.   Recommendations 

In  this  analysis,  the  assumption  used  is  the  elastic 
behavior  of  deformations.   As  stated  in  the  conclusion,  it  is 
still  believed  reasonable  to  use  the  elastic  theory,  but  for 
more  accurate  prediction,  local  plastic  deformation  phenomena 
should  be  included  in  the  analysis. 

If  the  elastic  assumption  is  to  be  used,  more  careful 

analysis  to  evaluate  the  value  of  6   and  C  are  recommended  in 

o 

future  research.   Also,  a  more  rational  approach  to  express 
C  is  recommended. 

In  future  research,  it  is  recommended  that  series  of 
experiments  be  conducted  to  measure  the  maximum  plate  defor- 
mation only,  and  express  C  and  8   in  terms  of  welding  and 
plate  dimension.   In  this  way,  the  computer  inputs  are  always 
available  for  any  prediction  of  deformation  for  the  two 
dimensional  panel  structures. 
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APPENDIX  B 
Example  of  Assembling 

As  an  example,  the  assembling  procedure  for  a  quarter 
part  of  the  plate  divided  by  two  finite  element  shown  in 
Figure  18a  will  be  discussed. 

For  the  element  1,  as  shown  in  Figure  18b ,  the  stiffness 

matrix  in  a  local  coordinate  can  be  written  from  equation 

(58a,  b) , 

k.   =  k  +  k  (A-l) 

~tl   ~   ~wx 

where  k   is  16  by  16  symmetric  matrix.   And  from 

1 
equation  (58d) ,  the  load  matrix  for  the  element  1  is, 

Q.   -  Qv  (A-2) 

where  Q,   is  16  by  1  matrix.   By  the  same  way  for  element  2, 

~tl 
the  stiffness  and  load  matrix  can  be  written  using  equation 

(58a,  b,  c)  and  (58d,  e) ,  respectively, 

k^   =  k  +  k    +  k  (A-3) 

~t_   —   ~wx   ~wy 

Q4.  =  Q  +  Q  (A-4) 

~t2   ~x   ~y 

where  k.   is  16  by  16  symmetric  matrix  and  Q    is  16  by  1 

^t2  ^2 

matrix. 

At  every  nodal  point  there  are  four  unknowns  (w,  w  ,  w  , 

x   y 

w   ),  therefore,  after  assembling  element  1  and  2,  the 
equations  have  to  be  solved  with  proper  boundary  and  symmetry 
conditions  are  24  by  24  simultaneous  equations,  which  is  in 
the  matrix  form: 
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6q  IK  q  -  Q]  =0 


(A-5) 


where  the  matrix  K  is  24  by  24  symmetric  matrix,  Q  is  24  by  1, 

T 
q  is  24  by  1,  and  q   is  1  by  24  matrix. 

All  the  quantities  in  equation  (A-5)  are  defined  by  the 

global  coordinate,  therefore,  the  assembling  procedure  is  to 

locate  the  quantities  defined  by  equation  (A-l,  A-2,  A-3,  A-4) 

to  the  proper  position  in  equation  (A-5) .   This  operation  can 

easily  be  established  comparing  each  nodal  points  numbered  in 

local  coordinates  to  that  of  global  coordinates: 

Element  1,  nodal  points       Element  2,  nodal  points 

Coordinates 

local  global 

1  3 

2  5 

3  .  6 

4  4 

And  noticing  that  numbering  of  the  generalized  displacement, 

q. ,  follows  the  relation, 

i  =  4m  -  g  (A- 6) 

j  =  4n  -  h  (A-7) 

where  i  and  j  represent  the  index  for  generalized  displacements, 

m  and  n  for  nodal  numbers,  g  and  h  are  set  to  be  3  for  w, 

2  for  w  ,  1  for  w  ,  and  0  for  w 

x         y  xy 

For  example,  one  of  the  stiffness  matrix  elements  at 

nodal  point  4  in  global  coordinate  can  be  calculated. 


Coordinates 

loca 

1          global 

1 

1 

2 

* 

3 

3 

4 

4 

2 

12 


>    X 


a.   Nodal  points  numbering  in  global  coordinate 


->* 


b.   Nodal  points  numbering  in  local  coordinate 


Figure  18   Elements  and  Nodal  Points  Numbering  in 

Global  and  Local  Coordinate  for  Sample  Assembling 
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From  equation  (A- 6)  and  (A- 7) , 

i  =  16  -  g  (A-8) 

j  =  16  -  h  (A- 9) 

But  the  same  node  4  becomes  node  3  at  element  1  and  node  4 

at  element  2  in  local  coordinate. 
For  element  1 , 

(i)1  =  12  -  g  (A-10) 

(j)1  =  12  -  h  (A-ll) 

and  for  element  2 , 

(i)2  =  16  -  g  (A-12) 

(j)2  =  16  -  h  (A-13) 

Therefore,  if  g  -   1  and  h  =  2,  then  from  equation  (A-8) , 
(A-9) ,  and  (A-5) , 

Kij  =  K15,14  (A"14) 

From  equations  (A-10) ,  (A-ll) ,  and  (A-l) 

k.  (A-15) 

1 

11,10 

and  from  equations  (A-12) ,  (A-13) ,  and  (A-3) , 

k.  (A-16) 

2 

Z  15,14 

Therefore,  element  of  the  stiffness  matrix  in  global 
coordinate  becomes , 

K15,14    =   \  -    f   kt,  (A"17) 

11,10  15,14 
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All  other  stiffness  matrix  element  can  be  expressed  by  the 

same  manner  described  above,  and  the  programs  of  assembling 

(13) 
the  stiffness  matrix  are  available  m  FEABL. 
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APPENDIX  C 
Description  of  Input  Data 

Input  data  for  the  programs  listed  in  Appendix  D  are: 

Cx'  Cy'  6o'  D'  tf  E'  Lx'  Lv'  M  and  N  as  described  in 
section  III.   In  addition  to  these,  control  card  for  calcula- 
tion of  different  cases  in  each  job  with  changing  above  input 
variables  are  considered.   By  using  this  calculation  control 
input  card  in  each  job,  series  of  calculation  with  changing 
some  of  the  input  variables  can  be  calculated,  and  the  defor- 
mation behavior  with  respect  to  these  changing  variables  can 
easily  be  visualized  by  plotting  the  deformation  versus  these 
variables . 

Detailed  input  data  card  arrangement  and  format  are 
listed  in  Table  4 . 


Table  4 
Input  Data  Layout  Form 
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Card    Input   Program 
Column  Format    Symbol 


Definition. and  Unit 


Sample 
Input 


Card  Number  1 


1-5 


15 


Card  Number  2 


1-15   E15.7 


NCASE 


LX 


16-30 

E15.7 

LY 

31-35 

15 

M 

36-40 

15 

N 

Number  of  cases  involved 
in  each  job 


Length  of  1/4  plate  in 
x,  inches 


14 


16.0E+00 


Length  of  1/4  plate  in  y   10.0E+00 
Number  of  elements  in  x        8 
Number  of  elements  in  y       5 


Card  Number  3 

1-12  E12.7       CX 

13-24  E12.7       CY 

25-36  E12.7  THETAZ 

37-48  E12.7  GNU 

49-60  E12.7  THK 

61-72  E12.7       E 


Equivalent  constrained 
welding  moment  as  in 
Figure  6,  lb-in/in  rad. 

Equivalent  constrained 
welding  moment  as  in 
Figure  6,  lb-in/in  rad. 


Free  joint  angular 
change,  8Q,  radian; 

Poisson's  ratio,  D 


Plate  thickness,  t, 
inches 

Young's  modulus,  E, 
lb/in.2 


4.20E+04 
4.20E+04 

55.0E-03 

0.30E+00 
0.375E+00 

30.0E+06 
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APPENDIX  D 
Listing  of  Programs 
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